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Abstract 

(N 

J> ■ We establish the existence of local, covariant time ordered products of local 

OP , Wick polynomials for a free scalar field in curved spacetime. Our time ordered 

products satisfy all of the hypotheses of our previous uniqueness theorem, so our 
construction essentially completes the analysis of the existence, uniqueness, and 
renormalizability of the perturbative expansion for nonlinear quantum field theories 
in curved spacetime. As a byproduct of our analysis, we derive a scaling expansion 
of the time ordered products about the total diagonal that expresses them as a sum 
O^' of products of polynomials in the curvature times Lorentz invariant distributions, 

^ , plus a remainder term of arbitrarily low scaling degree. 
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1 Introduction 

In order to give a perturbative definition of a nonlinear quantum field theory in a globally 
hyperbolic, curved spacetime, it is necessary to define Wick polynomials and their time 
ordered products for the corresponding linear (i.e., non-self-interacting) field. In the case 
of a scalar field, a construction of these quantities was given recently by Brunetti, Fre- 
denhagen and Kohler and by Brunetti and Fredenhagen 0. However, these authors 
did not impose a locality or covariance condition on the Wick polynomials or their time 
ordered products. In fact, the Wick polynomials were constructed in by means of 
a normal ordering prescription with respect to an arbitrarily chosen Hadamard vacuum 
state. The Wick polynomials defined in this manner thereby possess an undesireable non- 
local dependence upon the choice of this vacuum state. Since no locality or covariance 
condition was imposed on the construction of time ordered-products of Wick polynomials 
in P] — and, indeed, such conditions could not have been imposed since the Wick poly- 
nomials used in ||^ were not local, covariant fields — the renormalization ambiguities were 
found to involve coupling functions rather than coupling constants. 
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In a recent paper [Q, we introduced the notion of a local, covariant quantum fielcQ 
and we then imposed the requirement that the Wick polynomials and their time ordered 
products be local, covariant quantum fields. We also required that these quantities have 
a suitable continuous and analytic dependence upon the spacetime metric and have a 
suitable scaling behavior under scalings of the metric. (These latter notions are well 
defined only for local, covariant fields.) In addition, we required the Wick polynomials 
and their time ordered products to satisfy various additional properties, namely suitable 
commutation relations with the free field, a microlocal spectral condition, and (for the 
time ordered products) causal factorization and unitarity conditions. We refer the reader 
to for the precise statements of all of our conditions as well as a complete explanation 
of the algebraic framework within which our conditions were formulated. 

In [|T^ , uniqueness theorems were proven for both the Wick polynomials and their time 
ordered products. For the Wick polynomials, we showed that any two constructions that 
satisfy all of the above conditions can differ at most by a suitable sum of products of cur- 
vature terms of a definite scaling dimension multiplied by lower order Wick polynomials of 
the appropriate dimension. In particular, this implies that the ambiguity in defining Wick 
polynomials up to a given finite order is uniquely characterized by only a finite number 
of parameters. A similar uniqueness result was obtained for the time ordered products, 
thereby establishing that the ambiguity in defining these quantities up to any given finite 
order also is characterized by only a finite number of parameters. We then showed that 
Xip'^-theoTj in curved spacetime is renormalizable in the sense that the ambiguities arising 
in the perturbative definition of this theory correspond precisely to the (finite number of) 
parameters appearing in the classical Lagrangian (provided that the possible curvature 
couplings of the appropriate dimension are included in this Lagrangian). Again, we refer 
the reader to |]12| for the precise statements and proofs of these results. 

The above uniqueness theorems, of course, do not address the issue of whether there 
actually exists a construction of Wick polynomials and their time ordered products that 
satisfies all of our requirements. As already noted above, in Wick polynomials were 
constructed via a normal ordering prescription, but they fail to satisfy our requirement 
of being local and covariant. However, this deficiency can be repaired in a relatively 
straightforward manner by replacing the normal ordering prescription with respect to a 
(nonlocally defined) Hadamard vacuum state by a point-splitting prescription based upon 
a locally and covariantly defined Hadamard parametrix. It was proven in |[T^ that such 
a construction does indeed satisfy all of our requirements, thus establishing the existence 
of local, covariant Wick polynomials. 

One might hope that the construction of time ordered products given in could 
be similarly modified to yield local, covariant fields that satisfy all of our requirements. 
However, it is not at all obvious how to do this. As in (and as will be explicitly seen 
in subsection 3.1 below), the essential difficulty in defining time ordered products arises 

^ A general notion of local, covariant quantum fields has been given by M . 
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from the extension of certain multivariable distributions to the total diagonal; it is here 
that regularization/renormalization is needed. The usual momentum space methods of 
regularization are inapplicable in a curved, Lorentzian spacetime, but the Epstein-Glaser 
prescription |jlO| is well defined . However, this prescription involves the modification of 
test functions by the subtraction of their truncated Taylor series multiplied by a "cutoff 
function". The introduction of such a cutoff function makes the prescription inherently 
nonlocal. Consequently, the time ordered products defined by this prescription will fail 
to be local, covariant fields. 

A similar difficulty with the Epstein-Glaser prescription occurs in Minkowski space- 
time, where the introduction of the cutoff function makes the prescription fail to be 
Lorentz invariant. However, in Minkowski spacetime, a cohomology argument can then 
be used to establish existence of a satisfactory Lorentz invariant prescription |T^. We 
have not been able to generalize this argument to curved spacetime. For this reason, the 
issue of existence of local covariant time ordered products was left open in [|I^ . 



The main purpose of this paper is to prove the existence of local covariant time or- 
dered products, thereby essentially completing^ the perturbative construction of nonlinear 
quantum field theory in curved spacetime. The basic idea of our construction is as follows. 
As already indicated above, our task is to extend certain distributions on M"+^ \ A„+i 
to all of M""*"^ in a local, covariant manner, where M denotes the spacetime manifold, 
M"+i = x"+^M, and A„+i denotes the total diagonal of M"+\ 

A„+i = . . . ,x) I X G M}. (1) 

The key idea which enables us to accomplish this is to analyze the scaling behavior of the 
unextended distributions near the total diagonal. To do so, we first introduce n "relative 
coordinates" y and show that we can view each unextended distribution as a distribution 
in y for each fixed a; G M (i.e., our distribution in (n + 1) variables can be viewed as 
a distribution in the n relative coordinates that is parametrized by the point x on the 
total diagonal). We then show that near the total diagonal, each unextended distribution 
in question can be written as a finite sum of terms together with a "remainder term" 
with the following properties: (i) The terms in the finite sum are products of curvature 
terms in x times distributions, m, in y that correspond to Lorentz invariant distributions 
in Minkowski spacetime^. (ii) The remainder term has a sufficiently low scaling degree 
under scaling of y. The distributions, m, may then be extended to the total diagonal by 
Minkowski spacetime methods, whereas the remainder term can be extended to the total 



^ As pointed out in fl^ , when defining Wick polynomials involving derivatives of the field, it is natural 
to require the vanishing of any Wick product that contains a factor of the wave operator applied to the 
field. This requirement was not imposed in so the issue of existence of Wick polynomials that satisfy 
this additional condition remains open. We are currently investigating this issue; see alsi 

•^For the Feynman propagator and its powers, these terms would correspond to the m 
expressions given in pi , since each u can be given a momentum space representation. 
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diagonal by continuity. The resulting extended distributions can then be shown to provide 
a definition of local, covariant time ordered products that satisfy all of our requirements. 

The paper is organized as follows. In section 2, we review our requirements on the 
definition of time ordered products. These requirements are the ones previously given 
in except that we have replaced the continuity requirement of under smooth 
variations of the metric with a smoothness requirement. Further discussion of our new 
smoothness requirement is given in Appendix A. 

In section 3, we reduce the problem of constructing time ordered products to that of 
extending certain scalar distributions to the total diagonal. In subsection 3.1, we proceed 
inductively in the number, n, of variables, and reduce the problem to the extension of 
the time ordered products in n + 1 variables to the total diagonal. In section 3.2, we use 
a local, covariant version of the Wick expansion to express these time ordered products 
in + 1 variables as sums of local Wick products times "c- number" distributions, 
In subsection 3.3, we then translate our requirements on the definition of time ordered 
products into requirements on the extensions of the distributions, t", to the total diagonal. 

Section 4 is devoted to obtaining the desired extension of In subsection 4.1, we 
introduce "relative coordinates", y, and then derive our scaling expansion of with the 
properties indicated above. (Some properties of the distributions occurring in the scaling 
expansion are obtained in Appendix B.) The scaling expansion is then used to extend t° 
in subsection 4.2. Finally, in subsection 4.3, we show that the extended distributions, t, 
satisfy the properties listed in subsection 3.3, so that they define a notion of time ordered 
products satisfying all of the requirements of section 2. Some concluding remarks are 
given in section 5. 

We will restrict consideration here to the theory of a scalar field ip, but our basic meth- 



ods and results should be applicable to other fields. As in [0, for notational simplicity 
we restrict attention to time ordered products of Wick powers that do not contain deriva- 
tives if. However, our results should extend straightforwardly to time ordered products 
involving derivatives of the field, subject to the caveat mentioned in footnote ^ above. In 
addition, for notational simplicity in treating the scaling behaviour, we restrict consider- 
ation to the massless case, so that the free theory contains no dimensional parameters. 
Again, our results can be straightforwardly generalized to the case where dimensional 
parameters are present. 



Notation and Conventions. Our notation and conventions are the same as in [|T^ 
In particular, we define the Fourier transform on by u{k) = {27i)~^^'^ J u{x)e^^^^d^x. 
Multi-indices are denoted by a = (ai, . . . , a^) G N^. If a is an m-dimensional multi- 
index, then we also use standard notations such as |a| = ^a^, x° = x"^ . . . x^'" and 
8°" = g^."!^' g^Qm • We also use the "constant convention", meaning that we use the same 
symbol C for possibly different numerical constants in a chain of inequalities. The space of 
compactly supported smooth functions on a space X with values in the complex numbers 
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is denoted by 'P(X) and the space of smooth functions on X (not necessarily of compact 
support) by S{X). (For the definition of the topology on these spaces, see e.g. [|l^, Chap. 
V].) The corresponding topological dual spaces of distributions are denoted by V{X) 
respectively S'{X). The elements in S'{X) are the distributions of compact support. The 
wave front set ||14| of a distribution u is denoted by WF(m) and its analytic wave front 
set [jl4| (see also Appendix A) is denoted by WFa(m). 



2 Required Properties of the Time Ordered Products 

For the theory of a free scalar field, f{x), on an arbitrary globally hyperbolic spacetime, 
(M, g), we previously defined |T^ an "extended Wick-polynomial algebra", W(M, g), 
which generalizes the construction of Diitsch and Fredenhagen to curved spacetimes. 
This algebra is sufficiently large to contain elements corresponding to all Wick powers, 
(p^{x), (as distributions on compactly supported test functions on M) and their time 
ordered products 

T = T{^'^{x,)...^'"{xn)), (2) 



(as distributions on compactly supported smooth test function on M"). In ||T2[ we imposed 
a set of requirements on both and T that uniquely determined these quantities up to 
certain well specified renormalization ambiguities. In [|1^, we also constructed Wick 
products satsifying all of our conditions, so in this paper we will view these quantities as 
known. Our task here is to construct time ordered products of Wick powers that satisfy 
the following list of requirements, which — apart from the smoothness condition T4 (see 
remark (1) at the end of this section) — correspond to the requirements previously given 
in m 



Tl Locality/Covariance. The time ordered products are local, covariant fields, as 
defined in [||]. 



T2 Scaling The time ordered products scale "almost homogeneously" under rescalings 
g —>■ A^^g of the spacetime metric in the following sense. Let $ be a local, covariant 
field in n variables, and let be the rescaled local, covariant field given by iSA$[g] = 
A~'^'^(jA$[A~^g], where ax : >V(M, A~^g) —>■ >V(M, g) is the canonical isomorphism defined 
in 



12|. The scaling dimension, dq,, of a local covariant field is defined as 



= sup{(5 G R I lim X'^Sx^^ = 0}. (3) 
The scaling requirement on the time ordered product is then that 

N 

X-dT^^T = T + ^ In'^ A ^h, (4) 

h=l 
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where (It = Yl^i^ some natural number and where are local, covariant fields with 
scaling dimension (It which have fewer powers in the free field than T. 



T3 Microlocal Spectrum condition. Let ut be any continuous state on VV(M, g), 
so that, as shown in uj has smooth truncated n-point functions for n ^ 2 and a 
two-point function U2{x,y) = u{(p{x)(p{y)) of Hadamard from, i.e., WF(a;2) C C+(M, g), 
where 

C+(M, g) = {(xi, h; X2, -k2) e T*M^ \ {0} | (xi, h) ~ (x2, ^2); h e {V+),,}. (5) 

Here the notation (xi, ki) ~ (x2, ^2) means that Xi and X2 can be joined by a null-geodesic 
and that ki and k2 are cotangent and coparallel to that null-geodesic. (V""'")^: is the future 
lightcone at x. Furthermore, let 

n 

UT{Xu...,Xn)=Uj{T{[[^''^{x,))). (6) 



Then we require that 



WFK)cCT(M,g), (7) 



where the set Ct(M, g) C T*M" \ {0} is described as follows (we use the graphological 
notation introduced in @, ||): Let G{p) be a "decorated embedded graph" in (M, g). By 
this we mean an embedded graph C M whose vertices are points xi, . . . ,x„ G M and 
whose edges, e, are oriented null-geodesic curves. Each such null geodesic is equipped 
with a coparallel, cotangent covectorfield Pe- If e is an edge in G{p) connecting the points 
Xi and Xj with i < j, then s(e) = i is its source and t{e) = j its target. It is required that 
Pe is future/past directed if Xs(e) ^ J'^{xt(e))- With this notation, we define 

Ct{M, g) = |(xi, ki, . . Xn, kn) G T*M'^ \ {0} | 3 decorated graph G{p) with vertices 

xi, . . . , Xn such that ki = pe — Pe Vi|. (8) 

e:s{e)=i e:t(e)=i 

T4 Smoothness. The functional dependence of the time ordered products on the space- 
time metric, g, is such that if the metric is varied smoothly, then the time ordered products 
vary smoothly, in the following sense. Consider a smooth one parameter family of metrics 
g*^^\ let T^^^ be a corresponding family of time ordered products, and let C':^^ be given by 
eq. @ for this family of metrics. Furthermore, let o;^^-' be a family of Hadamard states 
with smooth truncated n-point functions (n 7^ 2) depending smoothly on s and with 
two-point functions ^2"^ depending smoothly on s in the sense that (see Appendix A) 

WF(cu2) C |(s, p; xi, X2, A;2) G r*(Rx M2)\ {0} I (xi,A;i;x2,A;2) GCj!^}, (9) 
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where the family of cones C^^ is defined by eq. (|D in terms of the family g'^'^^. Then we 
require that the family of distributions given by 

n 

Ms, Xu...Xn)= J'\T^'\l[ ^"'{Xi))) (10) 

i=l 

depends smoothly on s with respect to ■* in the sense that 

WF{ut) C p; xi, fci;... ; fc„) gT*(Mx Ar)\{0} I (xi,fci;...;x„, G 4'^}. 

(11) 

T5 Analyticity. Similarly, we require that, for an analytic one-parameter family of 
analytic metrics, the expectation value of the time ordered products in an analytic family 
of states varies analytically in the same sense as in T4, but with the smooth wave front 
set replaced by the analytic wave front set. 

T6 Symmetry. The time ordered products are symmetric under a permutation of the 
factors. 

T7 Unitarity. We have T* = T, where T is the "anti-time-ordered" product, defined 
as 

f (^^H^i) • • • = i-^r^'nU nw ^'^{x,)), (12) 

/lU---U/j={l,...,n} ieh i&Ij 

where the sum runs over all partitions of the set {1, . . . , n} into disjoint subsets Ii, . . . 

T8 Causal Factorization. In the case of a single factor, we require that T{ip^{x)) = 
ip^{x). For more than one factor, we require the time ordered product to satisfy the 
following causal factorization rule, which reflects the time-ordering of the factors. Consider 
a set of points (xi, . . . , Xn) G M" and a partition of {1, . . . , n} into two non-empty disjoint 
subsets I and /'^, with the property that no point x-i with i G / is in the past of any of 
the points Xj with j G that is, Xi ^ J~{xj) for all i G / and j G P. Then the time 
ordered products factorize in the following way: 

T = T{l[^''^{x,))T{l[^'^{x,)). (13) 
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T9 Commutator. The commutator of a time ordered product with a free field is given 
by lower order time ordered products times suitable commutator functions, namely 

n 

[T{^''{x{) . . . ¥.'"(x„)), ^{y)] = i k^Aix,, y)T{^'K^i) ■ ■ ■ . . . (14) 

1=1 

where A is the causal propagator (commutator function), defined as the difference between 
the advanced and retarded fundamental solutions of the Klein-Gordon equation. 



Remarks. (1) In our paper [12], we defined a notion of the continuous variation of a 



local covariant field under smooth variations of the metric, and we imposed this as a 
requirement on Wick powers and their time ordered products. We have replaced this 
requirement here with condition T4, which requires smooth (rather than continuous) 
dependence of the fields. It is easy to verify the the uniqueness results of [|l^ as well as 
the existence result of for Wick powers go through without any essential change if 
the continuity requirement imposed there is replaced by condition T4. We prefer to work 
with condition T4 here because it is a much simpler condition to state, it is more general, 
and it closely parallels the analyticity requirement T5 that was previously imposed in 
[0 . Further discussion and explanation of conditions T4 and T5 is given in Appendix A. 

(2) The microlocal spectrum condition is the same condition as formulated in ||^. It may 
be motivated by the fact that for noncoinciding points, u(T{Y[f^'{xi))) can be expressed 
in terms of Feynman graphs. A line in such a graph represents a Feynman propagator, 

ljf{x, y) == uo{T {ip{x)ip{y))) = uj2{x, y) — iA^'^^(x, y), whose wave front set off the diagonal 
is given by |1T8| 

W¥{ujf) = {(xi, fci; X2, -k2) I (xi, ki) ~ (x2, ^2); h e {y^)x, ^ X2 G J^(xi)}. (15) 



For non- coinciding points, the form of WF(u;t) follows from (^) and the rules for calculat- 
ing the wave front set of a product of several distributions, see e.g. [p!4| , Thm. 8.2.10]. For 
coinciding points, the form of WF(c(jj') reflects the usual energy momentum conservation 
rules. On the total diagonal, A„, the microlocal spectral condition reduces to 

^NY{UT) U„^T(A„,). (16) 

where the notation "±" means the following. If F C T*X with X a manifold and Y d X 
a smooth submanifold, then F fyX TY means that for any (y, k) & F \y and any 
(y, v) G TY we have that kaV°' = 0. 

(3) The "connected time ordered product" , T'^, of n Wick-monomials is deflned in terms 
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of the time ordered product by 

A" I 

/iU---U/j={l,...,n} iG/i ielj 

where the fi are test functions of compact support and S{f) is the formal S'-matrix for 
the Lagrangian ^{x) = J2i fii^)'^'''i^)j 

s{f) = J2-J nn^i) . . . ^M) /ig(xi) . . . (17) 

Our unitarity condition, T9, is equivalent to the condition T'^* = (— l)"+iT^ on the 
connected time ordered product. 

(4) For Minkowski spacetime, condition T9 was given in |P, IH, where it was shown to be 
equivalent to the familiar Wick-expansion of the time ordered products (see subsection 
3.2 below). 

Our task is to construct time ordered products of Wick powers that satisfy conditions 
T1-T9. We shall proceed inductively in the number of factors, n, appearing in the time 
ordered product (H). By condition T8, for n = 1 the time ordered products are just 



the Wick powers, which were already constructed in |12|. Therefore, we may inductively 
assume that time ordered products with properties T1-T9 have been defined for any 
number of factors < n. The goal is to construct from these the time ordered products 
with n + 1 factors. In the next section, we reduce the problem (in close parallel with the 
analysis of 0) to that of extending certain multivariable scalar distributions to the 
total diagonal. 



3 Reduction to the problem of extending certain scalar 
distributions to the total diagonal 

3.1 Construction of time ordered products up to the total diag- 
onal 

The key idea of causal perturbation theory is that the time ordered products with n + 1 
factors are already uniquely determined as algebra-valued distributions on the manifold 
^n+i ]2iinus its total diagonal A„_|_i by the causal factorization requirement T8 (see 
eq. (0)), once the time ordered products with less than or equal to n factors are given. 
Following 0, this can be seen as follows: 
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Let J be a proper subset of{l,2,...,n+l}, and let Ci be the subset of M""*"^ defined 

by 

Ci = {{xi,X2,...,Xn+i)\xi^ J^{xj) for alH e /, J G r}, (18) 

where P is the complement of /. It can be seen that the sets Cj are open and that the 
collection {Ci} of these sets covers the manifold M""*"^ \ A^+i. Let {//} be a partition of 
unity subordinate to this covering. On the manifold M"'^^ \ A„_|_i, we define the algebra- 
valued distributions T° by 

/C{l,...,n+l},7^0 

where 

Ti = n\{^\x,))T{\{^'^{x,)). (20) 

Using the causal factorization property T8 of the time ordered products with less or equal 
than n factors, it can be seen that the definition of T° does not depend on the choice 
of the partition {//}, so is well defined. Property T8 applied to the time ordered 
products with n + 1 factors then requires that the restriction of T to M^^^ \ A„+i must 
agree with T°. Thus, property T8 alone determines T up to the total diagonal, as we 
desired to show. 

We now claim that — assuming that time ordered products with less or equal than n 
factors have been defined so as to satisfy properties T1-T9 on M" — the fields T° with 
n + 1 factors automatically satisfy^ the restrictions of properties T1-T9 to M""*"^ \ A^+i. 
Condition T8 can be immediately seen to hold by virtue of the definition of T". The proof 
that properties Tl, T2, T6, T7 and T9 hold is relatively straightforward. A proof of the 
microlocal spectral condition, T3, can be given in exact parallel with reference 0. A 
generalization of this argument can be used to prove that the smoothness and analyticity 
conditions, T4 and T5, also hold. 

Our remaining task is to find an extension of each of the algebra-valued distributions 
T° in n + 1 factors from M"+^ \ A„+i to all of in such a way that properties T1-T9 

continue to hold for the extension. This step, of course, corresponds to renormalization. 
Condition T8 does not impose any additional conditions on the extension, so we need 
only satisfy T1-T7 and T9. However, it is not difficult to see that if an extension T is 
defined that satisfies T1-T5 and T9, then that extension can be modified, if necessary, 
so as to also satisfy the symmetry and unitarity conditions, T6 and T7. Namely, if the 

^Of course, if any T° failed to satisfy any of these properties on Af"+^ \ A„+i, we would have a proof 
that no definition of time ordered products could exist that satisfies T1-T9. 
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extension, T, of satisfied T1-T5 and T9 but failed to satisfy the symmetry condition, 
T6, we could define a new extension T' by symmetrization, 

^ ' Perm vr 

The so obtained extension of T° then satisfies T1-T6 and T9. Similarly, suppose the 
extension, satisfied T1-T6 and T9 but failed to satisfy the unitarity condition, T7, so 
that the corresponding connected time ordered product, T^, fails to satisfy T^* = (— 1)'^T^ 
(see remark (3) of section 2). Then we define T^' = |(T'^ + (— l)"T'^*) and redefine our 
extension by 

T' = T^\fl^''{x.))- Yl ^-^^Til[^''ix,))...Til[^'^{x,)). (22) 

/lU---U7j = {l,...,n+l} 

i>2 

This provides us with an extension of T° that satisfies T1-T7 and T9. 

Thus, we have reduced the problem of defining time ordered products to the problem 
of extending the distributions T° defined by (|19|) from M""*"^ \ A„+i to all of M"+^ so that 
properties T1-T5 and T9 continue to hold for the extension. In the next subsection, we 
will see that property T9 can be replaced by the requirement of a local Wick expansion 
for time ordered products. 

3.2 Reduction to a c-number problem via a local Wick expan- 
sion 

The next key simplification is to reduce the problem of defining the algebra valued distri- 
butions T to the problem of defining certain "c-number" distributions t. As in this is 
accomplished by means of a "Wick expansion" . The usual Wick expansion in Minkowski 
spacetime expresses time ordered products as a sum of normal ordered products with 
distributional coefficients. In the generalization to curved spacetime given in , the time 
ordered products are Wick expanded in terms of normal ordered products defined relative 
to an arbitrarily chosen quasi-free Hadamard state. However, such an expansion would 
not be useful here because the quasi-free Hadamard state — however it is chosen — has a 
nonlocal character. Consequently, the distributional coefficients occurring in the Wick 
expansion with respect to normal ordered products will fail to inherit the locality and 
covariance properties of the time ordered products themselves. 

For this reason, we will employ here a Wick expansion of the time ordered procducts 
with respect to the local, covariant Wick products : ip'^'^ . . . ip^" :h that were previously 



defined in as follows: Let H{x, y) denote the local Hadamard parametrix 

H{x,y) = Uix,y)a-' + V{x,y)\na, (23) 
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where f/, V are certain smooth functions defined in terms of the metric and the couphng 
parameters, cr is the signed squared geodesic distance and where the "iO" prescription for 
the singular terms is as for the two-point function in Minkowski space. (The formal power 
series defining V need not converge in smooth, non-analytic spacetimes, but a suitably 



modified convergent V can be used, as explained in |jT2|, Sec. 5.2].) Following |jT2|, Sec. 
5.2], we then define in some neighborhood, Un, of the total diagonal, A„, the algebra 
valued distributions 



(24) 



with \k\ = Yli^i- C)ur Wick-expansion is 

ni{v''{x,)) = J2 (%n..,Axu...,Xn) :/-^-Hxi)...V^'"-^"(x„):h, (25) 

where are c-number distributions on Un and where = Yl {f )- Note that 

the Wick-expansion formula (|25| ) is a different identity for different sets of exponents 
(fci, . . . , kn), but that the same coefficients tj^,,,j^ appear in each identity. 

We emphasize that since our local Wick-products ( |2^ ) are defined only on a sufficiently 
small neighborhood, Un-, of the total diagonal, our Wick-expansion will only make sense in 
this neighborhood. (This is in contrast with the Wick-expansion used in |^ which is based 
on a normal ordering prescription for Wick-products and therefore makes sense everywhere 
on M"'.) This fact, however, will not cause any complications for our constructions, since 
we will need the Wick-expansion only for the purpose of extending to the total diagonal. 

We claim now that any definition of time ordered products that satisfies requirements 
T3 and T9 must admit a Wick expansion of the form (P^D, with distributional coefficients 
satisfying 

WF(t,,..,,JcCT(M,g), (26) 

where Ct(M, g) is the set specified in (||). (Note that ( pG] ) implies in particular that the 
products of distributions implicit in our Wick-expansion formula actually exist and that 
the operator given by this formula defines — after smearing with a smooth test function — 
an element of our algebra W(M, g).) To prove this claim, we note that eqs. (pSj) and (^) 
hold trivially for the time ordered product T{ip) of a single free field. Let us now in- 
ductively assume that eqs. ( pS] ) and (^) have been demonstrated for all time ordered 
products of the form T{ip^^ . . . (p'^"), whenever \k\ = ki < d for some d > 1. We claim 
that they also hold for all multi-orders {ki, . . . , kn) with ^ki = d. To see this, we consider 
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the difference, 



1=1 



E 



jl---jn,\j\<\k\ 



t,,...,„(xi,...,x„) ly^'^^ ^Hxi) . . . y^'" '"ixn):H, (27) 



between the left side of eq. and the expression on the right side of that equation, but 
with the term tfci.,,A:„ 1 omitted in the sum. (Note that this is precisely the term in eq. (p5|) 
which is not already known by the induction hypothesis.) We now commute D with a free 
field (f. We use T9 to evaluate the commutator with the time ordered product and we use 
the similar commutation relation that holds for the local Wick products occurring in the 
sum. If this is carried out, then one finds that [D{xi, . . . ,Xn),f{y)] = 0. Since the only 
elements of our algebra W(M, g) that commute with all smeared field operators <f{f) are 
multiples of the identity |T^, Prop. 2.1], we thus find that D must in fact be given by a c- 
number distribution times the identity. We define tki...kn to be this c-number distribution. 
Now t = tki,„k„ can trivially be written as t = uj{D), for any Hadamard state, and each 
operator in the expression (p7|) for D satisfies f\ T3. Hence, condition T3 holds also for D, 
thus showing that WF(tfc^...fc^) C Ct(M, g). We have therefore completed the induction 
step, thereby establishing that the Wick-expansion holds for all multi-orders (fci, . . . , kn), 
provided only that T3 and T9 hold for the time ordered products. 

Conversely, if a definition of time ordered products has been given that admits a Wick 
expansion of the form (|25| ) with coefficients satisfying (^), then properties T3 and T9 
will hold as well in the neighborhood of the total diagonal on which the Wick expansion 
is defined. 

Since the distribution T° defined on M"+^ \ A„_,_i by (|T^) above satisfies properties 
T3 and T9 on M""^^ \ A„_|_i, it also admits a local Wick expansion of the form (p5|), i.e., 
on Un+i \ A„+i we have 

n+l 

Xi,...,Xn+l) ■.ip'''~^'{Xi) . . .ip''-+'-^-+'{Xn+l):H ■ 

(28) 



i = l jl-jn + l ^-^^ 



In the next subsection, we will reformulate the problem of extending the algebra-valued 
distributions to the algebra-valued distributions T in terms of the extension of the 
c-number distributions appearing in (|28|) to the c-number distributions t appearing in 



^For the terms in the sum in eq. ([27|), this follows from inductive hypothesis eq. ( Pq ) on the iji...j„ 
with |j| < together with the fact that uj{:ip''^ {xi) . . . tp''" (xn) -h) is smooth for all fci, . . . ,kn- 
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3.3 Reformulation in terms of the extension of t 

We return now to our inductive construction of time ordered products We assume that 
all time ordered products involving < n factors have been constructed so as to satisfy 
our assumptions T1-T9 and we consider an arbitrary time ordered product, T, in {n + 1) 
factors. As noted in subsection 3.1, property T8 will hold if and only if T is an extension 
to all of M"+^ of the distribution T° on M"+^ \ A„+i defined by (0). Since T° satisfies 
T1-T9 on M""*"^ \ A„+i, we need only check that our extension preserves these properties. 
As noted at the end of subsection 3.1, we actually need only check that T preserves 
properties T1-T5 and T9, since T8 does not provide any additional conditions on the 
extension and, by a suitable redefinition, it is straightforward to ensure that T6 and T7 
are satisfied. Furthermore, as shown in the previous subsection, we may replace property 
T9 by the local Wick expansion (^Sf). Thus, time ordered products satisfying all of our 
conditions will exist if and only if the c-number distributions t° on Un+i \ A„_|_i appearing 
in (1^) can be extended to distributions t on Un+i in such a way that the distribution T 
defined by ( |25| ) continues to satisfy properties T1-T5. It is straightforward to check that 
this will be the case if and only if the extensions t satisfy the following 5 corresponding 
conditions: 

tl Locality /Covariance. The distributions t are locally constructed from the metric 
in a covariant manner in the following sense. Let ip : N ^ M he a causality- preserving 
isometric embedding, and let / be a test function supported in a sufficiently small neigh- 
borhood in the total diagonal of N^~^^. Then we require that 

rt[gM]{f) = t[gN]{f), (29) 

where gM and gN are the metrics on M and respectively, so that 4'*gM = Sn- 

t2 Scaling. The distributions t = scale homogeneously up to logarithmic terms, 

in the sense that there is an A^ e N such that 

\-'t[X-h]=t[g] + J2^v,[g], (30) 

h=l 

where the Vh are certain local and covariant distributions, and where d = YlJi- 
t3 Microlocal Spectral Condition. WF(t) ^a^+i-L T(An+i). 

t4 Smoothness. Let g^*-* be a smooth family of metrics on M, depending smoothly 
on a parameter, and view t(s,xi, . . . ,x„+i) = )f:[g'-'*^](xi, . . . ,x„+i) as a distribution on 
M X Un+i- Then we require that 

WF(t) \rxA„+i'^ {{s,p;x,ki; . . .■,x,kn+i) I'^ki = 0, not all fci = 0}. (31) 
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t5 Analyticity. If g'^''^ is an analytic family of real analytic metric on Un+i, then item 
t4 holds with the smooth wave front set replaced by the analytic wave front set. 

Remarks. (1) If we apply the differential operator Xdx = a total of times to 
both sides of eq. pO|), then we obtain 

{Xdx-dft[\-'g]=0. (32) 

Moreover, if we apply A^a = only h < N times to both sides of ( pO] ) and set A equal 
to one afterwards, the we find that the local covariant distributions Vh are given by 

VH[g] = {Xdx-d)'^t[X-'g]\,=,. (33) 

In fact, equation (|3^) is actually equivalent to (|30D, as one can see by rewriting the 
differential operator Xdx — d as ■^'^g^-^"'^ and then integrating (|32D N times. 
(2) In formulating conditions t3-t5, we have taken advantage of the fact that on Un+i \ 
An+i, we have t = so t is already known to satisfy the wave front set conditions 
corresponding to T3-T5 on W„+i\A„+i. Consequently, we need only require t to satisfy the 
desired wave front set conditions on A^+i. Similarly, conditions tl and t2 also are already 
known to hold on Un+i \ A„+i, so we need only check that t satisfies these conditions in 
an arbitrarily small neighborhood of A„+i. 

In summary, in this section we have reduced the problem of defining time ordered 
products to the following question: Assume that time ordered products involving < n 
factors have been constructed so as to satisfy our requirements T1-T9. Define T° by (p!9|) 
and define the distributions on Un+i \ A„+i by (|28|). Can each be extended to a 
distribution t on Un+i so as to satisfy requirements tl-t5? 

4 Extension to the total diagonal 

Thus far, our analysis of time ordered products corresponds closely to that given in 
The primary difference in our assumptions is that we have imposed the requirement that 
time ordered products be local, covariant fields (see Tl) and that they satisfy certain 
additional requirements concerning scaling behavior (see T2), and smooth and analytic 
dependence on the metric (see T4 and T5). This has resulted in some important dif- 
ferences in our analysis as compared with In particular, as a consequence of the 
locality/covariance requirement, the Wick expansion of in terms of normal ordered 
products with respect to a quasifree Hadamard state is not useful, so instead we intro- 
duced a local, covariant Wick expansion in subsection 3.2. Nevertheless, all of the steps 
in the analysis given in section 3 above are in close parallel with the analysis of [Q. 
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As described at the end of section 3, our analysis will be completed if we can ex- 
tend the distributions t° to the total diagonal so that they satisfy properties tl-t5. As 
is well known from quantum field theory in Minkowski spacetime, straightforward at- 
tempts to extend to the total diagonal give rise to formal expressions that do not make 
sense as distributions. Therefore, one normally proceeds by introducing some means of 
"regularizing" these formal expressions and then extracting a well defined "finite part" 
(up to renormalization ambiguities). In Minkowski spacetime, most approaches to regu- 
larization/renormalization involve the use of Euclideanization and/or momentum space 
methods, neither of which have a natural generalization to (non-static) curved Lorentzian 
spacetimes. For this reason, the authors of ||^ employed the regularization procedure of 
Epstein and Glaser, which is "local" in the sense that it uses coordinate space methods 
that can be defined in a local region. 

Nevertheless, the Epstein- Glaser method is not local in a strong enough sense for our 
purposes, since we need to ensure that the renormalized time ordered products will be 
local, covariant fields. A key step in the Epstein-Glaser regularization procedure is the 
introduction of certain "cutoff functions" of compact support in the "relative coordinates" 
that equal 1 in a neighborhood of the total diagonal. Since the prescription for the 
extension of depends upon the spacetime geometry throughout the region where the 
cutoff functions are non-zero, the extension, t, at a point p G A„+i will not depend only on 
the metric in an arbitrary small neighborhood of p and, thus, will not depend locally and 
covariantly on the metric in the sense required by condition tl. There does not appear to 
be any straightforward way of modifying the Epstein-Glaser regularization procedure so 
that the resulting extension, t, will satisfy property tl. In particular, serious convergence 
difficulties arise if one attempts to shrink the support of the cutoff functions to the total 
diagonal. In addition, the cohomological argument of Jl^ also does not appear to admit 
a straightforward generalization to curved spacetime. 

Consequently, we shall proceed by a different route here. Our approach to extend t° 
to the total diagonal is motivated by the idea (essentially the "equivalence principle") 
that on sufficiently small scales a curved space "looks flat", and that the divergences 
of i° in curved spacetime should be of the same nature as the corresponding t° in fiat 
spacetime. However, this idea is not correct as just stated because a curved space is 
not actually fiat (no matter on how small a scale one looks). Although it is true that 
the leading order divergences of t° will be essentially the same as in fiat spacetime, in 
general there will be sub-leading-order divergences that are sensitive to the presence of 
curvature and are different from the divergences occurring for the corresponding in 
flat spacetime. Nevertheless, we will see in subsection 4.1 below that any local, covariant 
distributions that satisfies our scaling, smoothness, and analyticity conditions admits a 
"scaling expansion" about the total diagonal. This expansion expresses t° as a finite sum 
of terms plus a remainder term with the properties that (i) each term in the finite sum 
is a product of a curvature term times a distribution in the relative coordinates that 
corresponds to a Lorentz invariant distribution in Minkowski spacetime (which can be 
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extended to the total diagonal by Minkowski spacetime methods) and (ii) the remainder 
term admits a unique, natural extension to the diagonal by continuity. We shall thereby 
obtain an extension of in subsection 4.2. In subsection 4.3 we will then show that the 
resulting extension satisfies all of the required properties tl-t5. 



4.1 The scaling expansion 

As indicated above, the key step that will enable us to extend t° is to perform a scaling 
expansion of it about the total diagonal. However, a priori it is not even clear what this 
means, since t° is a distribution in n + 1 variables, and it is not clear what it would mean 
to perform any kind of "expansion" of a distribution about the 4-dimensional submanifold 

The first step in obtaining the scaling expansion for f is to show that it is possible to 
fix one of its n + 1 variables at a value x, and view it as a distribution in the remaining n 
variables, y, which play the same role as "relative coordinates" in Minkowski spacetime. 
In other words, writing 

x = xi, y = (xa, . . . (34) 

we show that the (unextended) distribution t° possesses a well-defined restriction to the 
submanifold 



a = Mx(f/"\(a;,...x)), (35) 

where [/ is a convex normal neighborhood of the point x G M. In Minkowski spacetime, 
this result would follow as an immediate consequence of translation invariance. In our 
context, this result follows from the microlocal spectral condition: Since property T3 
is known to hold for T*^, it follows that the wave front set of t° is contained in the set 
Ct- As can be seen from "energy- momentum conservation constraint" in eq. (|^), Ct does 
not contain any elements of the form (x, /c; ?/, 0). Since the conormal bundle, N*Cx of 
the submanifold is spanned precisely by such covectors, WF(t°) does not have any 
elements in common with N*Cx- That the restriction exists is thus ensured by Thm. 
8.2.4]. This restriction may be identified with a distribution on f/" \ (x, . . . , x). 

We now shall obtain our scaling expansion of . The basic idea is to expand the 
t''[g](a;, ■ ) at a fixed point x in terms of the metric and its derivatives at x. The individual 
terms in the so-obtained series will be seen to be given by sums of products of local 
curvature terms at x times Lorentz invariant distributions of the relative coordinates. 
The remainder for the suitably truncated series will not have this simple form, but will 
turn out to be regular enough to allow a unique extension. 

To begin, we choose a convex normal neighborhood U d M oi x and introduce Rie- 
mannian normal coordinates with respect to the metric g around x. These coordinates are 
constructed by using the exponential map to identify U with a subset of the tangent space 
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TxM, at X and then identifying T^M with Minkowski spacetime, M^, by an isomorphism 
e : T^M R^. Thus, the Riemannian normal coordinates of a point ^ & U are given by 

a,(0 = eo(expJ-l(OGM^ (36) 

However, when it is not hkely to cause confusion, we will slightly abuse the notation 
by denoting all quantities — i.e., the point, its Riemannian normal coordinates, and the 
corresponding vector in Minkowski spacetime — simply by ^. Similarly, the Riemannian 
normal coordinates of y = (x2, . . . ,Xn+i) (see eq. ( ^If ) above) will be denoted ax{y), but 
when it is not likely to cause confusion, we also will use y to denote the Riemannian 
normal coordinates of these points or the corresponding vector in R'^". 

The choice of isomorphism e : T^-M — > R'^ is equivalent to a choice of an orthonormal 
tetrad, e^, at the point x. Since any other orthonormal tetrad is of the form A'j^e'^ for 
some Lorentz transformation A, the Riemannian normal coordinates, ^, of a given point 
corresponding to the transformed tetrad at x are then given in terms of the original normal 
coordinates by A,^. Similarly, the Riemannian normal coordinates, y, of a point in f/" are 
obtained by Lorentz transforming the coordinates of each point individually by A, the 
result of which we shall denote as Ay. 

Now let g*^*-* be the smooth 1-parameter family of metrics on U whose coordinate 
components in Riemannian normal coordinates around x are given by 

9l:K0=9,uis0- (37) 

Note that if Xs is the map from U into itself given by ^ — > in Riemannian normal 
coordinates about x, then this family of metrics can be alternatively written as 

g^^^ = s-\:g. (38) 

Note also that the definition of the above family of metrics does not depend on any 
additional data besides the specification of the point x and the metric itself. In particular, 
it does not depend on our choice of tetrad at x. 

By a slight generalization of the microfocal argument given at the beginning of this 
subsection, it follows from the fact that T° satisfies properties T3 and T4 on U^~^^ \ 
An+i that t°[g*^^^](a;, ■ ) makes sense as a family of distributions on f/" \ (x, . . . , x) that 
is parametrized by {s,x). Furthermore, when smeared with a test function, /, in y, it 
follows that t^['g^^'^]{x, f) is smooth in {s,x). In addition, since differentiation does not 
increase the size of the wave front set, derivatives of with respect to s also make sense 
as distributions on f/" \ (x, . . . , x) that are parametrized by {s, x). Hence, for any k and 
any arbitrary, but fixed x G M, we can define a distribution on [/" \{x, . . . ,x) by 

(39) 

s=0 



r![g]{x,.) = ^^f [gW] (X,-) 
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It follows that for any given natural number m > 0, we have the following Taylor expansion 
with remainder: 

^° = E^-° + ^- (40) 

fc=0 

where 

1 fjm+l 

rMix, ■) = ^la- ^rj-^i' W'^] (^^ ■ ) ds. (41) 

Formula (|40|) is actually our desired scaling expansion of However, as it stands, (^) 
is merely an identity that would hold for any distribution in the variables {s,x,y) that 
satisfies suitable wave front set conditions. The important properties of this formula for 
our distributions are stated in the following theorem. 

Theorem 4.1. (i) t°(x, ■ ) and r^{x, ■ ) are distributions on f/" \ (x, . . . , x) which are 
locally constructed from the metric in a covariant way in the sense that eq. (^) 
holds for all diffeomorphisms that leave the point x invariant. 

(ii) We have the decomposition 

r',ix,y) = '£C{x)-a*y{y), (42) 

where the sum is finite. Here, C = C^^'"^^ denote the coordinate components 
of certain curvature tensors in Riemannian normal coordinates about x and the 
= n?, ,, are Lorentz-invariant tensor valued distributions defined on M^" with 
the origin removed, that is, 

<....,(a-) = a;:;---a>L.(-) (43) 

for any Lorentz-transformation A. The local curvature tensors C arise as a sum of 
monomials in Qab, Rated, ■ ■ ■ , V(ei . . . Vek_2)Rabcd- In the case considered here with no 
dimensional parameters, each monomial contains precisely k coordinate derivatives 
of the metric. 

(iii) r° and scale almost homogeneously under rescalings of the metric with degree d. 

(iv) The distributions scale almost homogeneously with degree d — k under coordinate 
rescalings in the sense that there exists an G N such that 

Slku' = 0, (44) 

where ^^^ = (Ee5Mf + P)^- 

(v) The scaling degree of r^(x, ■ ) is less or equal than d — m — 1, i.e., the distributions 
js^d-m-i+<5^o^^^^ ^ . -j^ viewed as distributions on W^^ \ via the pull-back by a^, tend 
to the zero distribution as A \ for all x and all 5 > 0. 
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Remark. We note that (iv) means that scales homogeneously with degree d — k 
under a rescaling of the coordinates, up to logarithmic terms. Namely, simple integration 
of (HP gives that 

X'^-'u^X ■ ) = ■ ) + E V ■ ^^^^ 

h=i 

This implies in particular that the scaling degree of vP at the origin is d — k. 

Proof. Item (i) follows directly from the fact that t° is local and covariant on its domain 
of definition. 

To prove (ii), we consider, first, the case where all the components of g in our Rie- 
mannian normal coordinates are polynomials in the Riemannian normal coordinates ^ in 
a neighborhood of x. Then we may characterize g by its components g^y at x together 
with the components of the coordinate derivatives, g^u,cTia2...i at x, only finitely many of 
which are nonzero. We may thus view t° as being a function of these quantities, and we 
express this by writing 

^°[g](a;, ■) = t^[g^,y, . . . , g^,u,a,a2...ai, ■ ■ ^{x, ■), (46) 

Now smear with a test function, /, on t/" \ (x, . . . , x). Since t°[g](x, /) depends smoothly 
on the metric — and hence is a smooth function of the finite number of variables 

on which it depends — we obtain 

rl[z]{xj) = [g(^)] (x, /) 



s=0 



^ [dfii^i ■ ■ ■ 1 S gf_ii/,ai(T2...cri ) • • • ] {x, f) 



=0 



ll+2l2+...mlm=k 



where [. . . ] stands for [g^u, 0, 0, . . . ]. 
We may rewrite this equation as 



T',{x,y) = J2C{x)-a:u\y), (48) 
where the C = C*^^" '*' are monomials in g^ij,aia2...amy which have the property that the total 



number of derivatives of (^^i, appearing in each C is equal to k, and where m° 
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are tensor- valued distributions on R^" minus the origin^, which are independent of g. 
Since we are working in Riemannian normal coordinates, the m-th coordinate derivates 
9fiu,aia2...iJm of the metric tensor at x can be rewritten as the coordinate components 
at X of a local curvature term that is polynomially constructed from the metric, the 
curvature tensor and its derivatives at x. Moreover, such a curvature term must involve 
precisely m derivatives of the metric. Hence, by our formula p7|), we conclude that 
(jiJ-i-'-fj-i corresponds to a curvature term C"'^"'"''- which arises as a sum of monomials in 
Rahcdi ■ ■ ■ , V(ei . . . ^ ek-2)Rabcd, Gach of which coutaius precisely k derivatives of the metric. 

We would next like to show that the distributions are Lorentz invariant. For this, 
we consider the diffeomorphism ip^ on U C M given by ^ ^ where A is a Lorentz 
transformation, and where the point ^ G f/ has been identified with its Riemannian normal 
coordinates about x. By definition, this diffeomorphism will leave the point x invariant, 
so we may apply item (i) to this diffeomorphism. From this, we get that 

J2 C^'-" (^X...,. (A • ) = E (^)^m • • • K<-u. ( ■ )• (49) 

Since this holds for all metrics, this means that must be Lorentz invariant. This proves 
(ii) for all metrics whose components in Riemannian normal coordinates are polynomials 
in the Riemannian normal coordinates ^ in a neighborhood of x. 

Now consider an arbitrary smooth metric g. In a compact neighborhood, K, of x, we 
can, for each n, find a metric q^"-* that is polynomial in ^ and is such that everywhere 
within K we have \gfj.u,aia2...am ~ (llj^,(Tia2...a,n\ < 2""' for all m < n. Let ip : R ^ [0,1] 
be a smooth function with support in [—1, 1] satisfying '?/'(— f) = "ipiv) and also satisfying 
1 — "ipiv) = "ipi^ — v) for all V G [0, 1]. Set h*^°^ = g and for s ^ but in a sufficiently small 
neighborhood of define 

h(^) = ^^(|l/s| -n)q('^) (50) 

n 

(Note that at each s, there can be at most two terms in this sum which are nonvanishing.) 
Then it is straightforward to show that h'-*-' is a one parameter family of smooth metrics 
that depends smoothly on s. Consequently, each r°[h(*)](x, ■) varies smoothly with s. 
However, we have already proven that eq. (^2]) holds for all s 7^ 0. By the smoothness 
property t3 applied to it follows that eq. (|^) continues to hold at s = 0, thus proving 
property (ii) for an arbitrary smooth metric g. 

Property (iii) is a direct consequence of the fact that t° satisfies the scaling property 
t2. To see this, we note that 

(Xd, - dfr'.lX-'g] = d', {\d, - dff [A~2g(^)] 1^^^ = 0, (51) 

^Actually, in the above construction the distributions u'^ are automatically defined only in the neigh- 
borhood of the origin in R**" (minus the origin itself) corresponding to the neighborhood, U C M on 
which the Riemannian normal coordinates are defined. However, by modifying g outside of a neighbor- 
hood of the origin if necessary, we may assume without loss of generality that the Riemannian normal 
coordinates are globally defined and that u° is defined everywhere on R''" \ 0. 
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since t° satisfies t2. This establishes (iii) for r^. That satisfies (iii) then follows 
immediately from eq. (0). 

To prove (iv), we note that (i) implies that 



= xUd'^t' [g(^^)] = X'xUr'.ig] (52) 



By eq. (|5TD , the differential operator {Xdx — d)'^ annihilates the left side of eq. ([5^). This 
implies that 

= {Xd, - df X'xUr'kls] = X' {Xd^ -d + kf xUr![s]. (53) 
Substituting the decomposition of r° into the expression on the right side, we obtain 

= ^C(x) ■(A9A-t^ + A;)^M°(A-) 

= J2C{x).Sl,u\X-). (54) 

Since this holds for arbitrary metrics g, it follows that that S^_}^u^ = 0, as we desired to 
show. 

In order to establish the estimate on the scaling degree for r^, item (v), we first use 
eq. d^ ) along with the same arguments as in (iv) to write 

rl{x,X-) = ^ / il-f.rdrh'[X'g'^^^]{x,-)\ df.. (55) 
Using the fact that t° satisfies property t2, we have (see eqs. (^) and (|33D ) 

N-l 

dx, A ■ ) = A"+i-^ J2 ^ ^°(^' ^' ■ ), (56) 

1=0 

where 

i^KKx, ■ ) = ^ - f^rdr\vd. - dYt' [v'g^^^] (x, ■ )L=,,,.=i^/^- (5T) 

If / is a smooth test function on f/"" whose support does not contain the point {x, . . . ,x), 
then by wave front set arguments similar to those given above, it follows from the fact 
that T° satisfies conditions T3 and T4 that the quantities ip'j^lXyX, f) are smooth in A in 
a neighborhood of zero. This immmediately implies (v). □ 
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Remarks. (1) As stated in property (v) of the above theorem, if we carry the scahng 
expansion, eq. (^O]) , to higher order (i.e., larger m), the remainder term rj^ will have 



a lower scaling degree. However, it should be noted that the wave front set of t is 
determined by the null geodesies of the curved spacetime metric g whereas the wave front 
set of each is similarly determined by the null geodesies of the fiat spacetime metric 
associated with the exponential map at x. Since the null geodesies of these two metrics 
do not, in general, coincide (with the exception of the null geodesies passing through x 
itself), it is clear that remains fundamentally distributional in nature no matter how 



large m is chosen. It also should be noted that it is not claimed in Thm. |4.1| that r„ 
converges to zero in any sense (even for an analytic spacetime) as m — >■ oo. Thus, eq. (^ 
should be viewed only as a "scaling expansion" with the properties specified in Thm 
not as a convergent power series. 

(2) If we combine eqs. ( ^0]) and (^2]), we obtain an expansion of of the general form 



t°(a;, y) = J2 C'(^) " «>°(z/) + r^ix, y). (5^ 



If the terms in the sum in (|58| ) are ordered by the engineering dimension of the curvature 
terms, C, then the first term in the expansion has C = 1 and the corresponding distri- 
bution is the "scaling limit" at x of the distributions t° in the sense of Fredenhagen 



and Haag |]rT|. The higher order terms in the expansion then give corrections to the the 
scaling limit, organized in powers of the curvature tensor and its derivatives. If dimen- 
sionful parameters are present in the theory, then the scaling expansion will be organized 
in terms of products of powers of the curvature and the dimensionful parameters. Our 
scaling expansion is also closely related to the "momentum space representation" of the 
Feynman propagator and its powers (see remark (3) below) given in |^ , since the Lorentz 
invariant distributions, vP, on Minkowski spacetime occurring in our expansion can be 
given a momentum space representation. 

(3) For the Feynman propagator and its powers, the scaling expansion can be explicitly 
calculated from known properties of the Hadamard expansion. We will illustrate this 
with two examples. The first example is the simplest nontrivial time ordered product, 
r°((^(a;) </?(?/)). Its Wick-expansion is given by 

T\^{x)^{y)) = ■.^{x)^{y):H +Hp{x,y)l, (59) 

where Hp = H — iA'^'^^ is the "local Feynman parametrix", where H is the Hadamard 
parametrix, eq. (|23D, and A'^'^^ is the advanced Green's function. Thus, the only nontrivial 



distribution t° occurring in this expansion is 

t°(x, y) = Hp{x, y) = U{x, y){a + 10)"^ + V{x, y) \n{a + iO), (60) 

where U and V are as in the Hadamard parametrix (see eq. (^)). The first few terms, 
Tj^, in the scaling expansion for = Hp are easily found from the expansions for U and 
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V given in ^ and many other references. Modulo an overall constant, one finds 



r^ix,y) = (r^^.ee'^ + iO)-i 
T^{x,y) = 

where, as above, ^'^ denotes the Riemannian normal coordinates of y relative to x. Thus, in 
this example, our scaling expansion corresponds to the usual short distance approximation 
to the singular part of the Feynman propagator (see, e.g., @). 

Our second example is the time ordered product T^{(p'^{x){p'^{y)). Its Wick-expansion 
is given by 

T\^\x)ip\y)) = ■.^\x)v\y):H +2Hp{x,y) ■.^{xMy):H +HF{x,y)H. (61) 

The only new arising in this expansion is the "fish graph", t° = Hp, a solution to the 
renormalization of which was found by B. S. Kay ||l5l prior to the commencement of the 
present work and played a role in the development of the present work. As can be seen 

r° 



from the above expansion for Hp, the first few coefficients, r^, for the fish graph are. 



modulo an overall constant 



r?{x,y) = 

r°(x, y) = lR'"'{x)Up{v,ueC + iO)-' - ^R{x){v,ueC + iO)"' Hv.pC^ + iO). 

It is easily seen that in both examples, the distributions are local, covariant distri- 
butions of the form claimed in (ii) — i.e., they are sums of terms of the form C{x) ■ a*u^{y) 
with a Lorentz-invariant Minkowski space distribution — and satisfy the scaling prop- 



erties specified in Thm. 



(4) The above scaling expansion was carried out for the scalar distributions t^. It is 
straightforward to check that it also holds for the extended distributions t that will be 
defined in the next subsection. Much more generally, it should be possible to perform a 
similar scaling expansion for arbitrary local covariant fields that satisfy appropriate wave 
front set properties. This should yield a generalized operator product expansion in curved 
spacetime. We are currently investigating the properties of such an expansion. 

4.2 Extension of t^[g] 

Theorem 4.1 of the previous subsection provides the necessary machinery to achieve our 
goal of extending t° in such a way that properties tl-t5 are satisfied. The basic idea is 



simply to suitably extend each term in the scaling expansion, eq. (^Ol) . Each in that 



equation is of the form (H3) and hence can be extended to the total diagonal by extending 
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the Minkowski spacetime distributions to the origin. This can be achieved by standard 
methods used in Minkowski spacetime. On the other hand, if m is chosen sufficiently large, 
the remainder term will have sufficiently low scaling degree that it can be extended 
to the total diagonal by continuity. The proof that the so-obtained extension t satisfies 
properties tl-t5 will be given in the next subsection. 

The key result needed to extend each is the following: 

Lemma 4.1. Let = fniv) with y = (^i,---,^n) be a Lorentz invariant tensor- 
valued distribution on W^^ \ which scales almost homogeneously with degree p under 
coordinate rescalings, i.e., 

M° = for some natural number A^. (62) 

where 5*^ = (^4^c?/c?^f + p)^ ■ Then m° has a Lorentz invariant extension, u, to a 
distribution on M^" which also scales almost homogeneously with degree p under rescalings 
of the coordinates. 

Proof. We will first extend using the Epstein-Glaser prescription. This extension need 
not satisfy either the scaling or Lorentz invariance properties. However, we will show 
that the extension can be modified, if necessary, so as to scale almost homogeneously^] 
with degree p. We will then show that the resulting extension can be further modified, if 
necessary, so as to be Lorentz invariant while retaining the almost homogeneous scaling 
with degree p. 

Choose an arbitrary smooth function w of compact support on W^"^ which is equal to 
one in a neighborhood of the origin. For any test function / G P(]R^") we set 

{Wf){y) = f{y)~w{y) y''dj{0)/a\, (63) 

\ct\<p—An 

where we use the usual multi-index notation. It follows from = that has scaling 
degree p, so by Thm. 5.3], we can define an extension, u, of to M'^" by setting 

uif)=u\Wf). (64) 

It follows that the scaling degree of m is p 0, Thm. 5.3], but it need not hold that u scales 
almost homogeneously with degree p, i.e., there is no guarantee that Sp^u = for some 
natural number M. However, one can calculate that 

WS^f{y)-S^Wf{y)= ^ r{y)dj{0) (65) 

\a\<p—4n 

^For distributions with an exactly homogeneous scahng, this result has previously been obtained in 
p^ , Thms. 3.2.3 and 3.2.4]. Thus, our theorem generalizes this result to the case of almost homogeneous 
scaling. 
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for some smooth functions '?/'° whose support does not contain the origin. From this it 
follows immediately that 

\Q\<p—4n 

where c" = 

We now define a modified distribution u' by 



u' = u 



El -A r-Ki^(^a5. (67) 

|«|<p-4n-l ^'^ ' 



Using the fact that SpdaS = [p — An — \a\)da5, we find 

S^u'= J2 c"^-^- (68) 

\a\=p—4n 

If we apply the operator Sp to both sides of the above equation, then we get that Sp^^u' = 
0, because 

SpdaS = for |a| = p - 4n. (69) 

This means that u' is an extension of m° with the desired almost homogeneous scaling. For 
notational simplicity, we will drop the "prime" in the following and denote this modified 
extension as u. 

We now investigate the Lorentz transformation properties of u. Restoring the tensor 
indices on u, we find by a calculation similar to eq. (p^ ) above that for any test function 
/ G P(]R^"') and any Lorentz transformation. A, we have 

^......,(/)-a-...a;:>.,....x^(a)/)= Y1 ^m.....(W(/), (70) 

|a|<p— 4n 

where {R{A)f){y) = f{Ay) and the fni-^) complex constants, which would vanish 
if and only if the distribution u were Lorentz invariant. We now apply the differential 
operator S^^^ to both sides of the above equation. Since Sp is itself a Lorentz invariant 
operator, we have R{A)Sp = SpR{A). Therefore, since Sp~^^u = 0, the operator S*^"*"^ 
annihilates the left side of eq. ([70|), so we obtain 

\a\<p—4n \a\<p—4n 

It follows immediately that &^^...^j(A) = 0, except possibly when \a\ = p — An. Thus, we 
have 

u,,...pAf) - K\ ■ ■ ■ A>.,...,(i?(A)/) = 6;:-;,;7-(A)9., . ..d,^_,J{f) (72) 
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for all / and all Lorentz-transformations A. Using this equation, one finds the following 
transformation property for 6(A), 



6(AiA2) = 6(Ai) + D(Ai)KA2), (73) 

where we have now dropped the tensor-indices and where D denotes the tensor represen- 
tation of the Lorentz-group on (®']R^)* ® (®^~^'^]R^). It then follows by the cohomological 
argument given in |T^ that this relation implies that h can be written in the form 



6(A) = a - D{K)a VA, (74) 

where a is an element in (®'ffi^)* ® (®''~^"M^), not depending on A. This enables us to 
define the modified extension 

^W-,. = - «MV.:.;r"^-i • • • 9u,.,j, (75) 

where we have now restored the tensor indices. It is easily checked that u' is Lorentz in- 
variant and satisfies S^^^u' = 0. We have therefore accomplished the goal of constructing 
the desired extension of u^. 

□ 

Some analyticity properties of u and its Fourier transform that follow from its scaling 
behaviour are established in Appendix B. These results, however, will not be needed in 
our present analysis. 

We now can give our prescription for extending t". Let d denote the scaling degree 
of let m = c? — An, and consider the expansion eq. (^Ol). By theorem 4.1, each r° 
appearing in this expansion takes the form 



r',ix,y) = }_^Cix)-ay{y) (76) 



where the sum is finite. We extend r° to a distribution Tk on f/" by choosing an extension, 
u, of each that satisfies the properties of Lemma ^4.1| and defining 



nix, y) = Y^ C{x) ■ alu{y). (77) 

Although Tfc has been constructed as a distribution in y that is parametrized by x, it is 
straightforward to check that may also be viewed as a distribution jointly in x and y. 

On the other hand, we know by property (v) of Theorem 4.1 that the scaling degree 
of rjjj is less or equal to 4n — 1. Therefore we can apply Thm. 5.2] to conclude that 
r^{x, ■ ) has a unique extension, r^ix, ■ ) to all of f/" with the same scaling degree for 
any given point x. This extention is given by 

rUf) = lim r^(t9(^V), (78) 
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where d^^'' is a sequence of smooth functions with support in f/"+i \ A„+i, which are 
identically one outside neighborhoods U^li of A„+i, with U^li shrinking to A„+i as j 
goes to infinity. By the scaling properties of r^, this limit exists in the weak sense, and 
is independent of the particular choice of cutoff functions i!)^^^ (see Thm. 5.2]). Again, 
it can be shown that this extension defines a distribution jointly in x and y. 
Our extension, t, is then defined by 

m ^ 
fe=0 

Our remaining task is to show that t satisfies properties tl-t5. 



4.3 Proof that t satisfies properties tl— 15 

As we now shall show, it is relatively straightforward to prove that the extension, t, of t° 
defined by eq. (|79D above satisfies properties tl and t2. 

To show that tl holds, we note that the prescription for extending clearly is local in 
the appropriate sense. However, it is not immediately obvious that the prescription yields 
a covariant extension in the sense required by tl since the prescription involves ax, 
whose definition requires, in addition to the metric, a choice of a tetrad at x. However, 
since any other tetrad at x is related by a Lorentz-transformation, it follows immediately 
from the Lorentz invariance of the extensions u in (|77| ) that different choices of tetrad 
lead to the same distribution r^. It follows that each is locally constructed from the 
metric in a covariant way in the sense required by tl. 

In order to see that is local and covariant in the sense of tl, it is sufficient to show 
that rm['?/'*g] is equal to '?/'*rm[g] for any diffeomorphism ip on U. We already know that 
this is true off the total diagonal A,„_|_i, as the unextended distribution r° has this property. 
Thus, the difference between the two expressions must be a distribution supported on the 
total diagonal. Moreover, the scaling degree of this distribution must be less than 4n — l, 
by our choice m = d — 4n. It is well known that there are no such distributions apart 
from the zero distribution (essentially because the delta function and its derivatives have 
scaling degree > An). Therefore the difference must in fact be zero, showing that 
satisfies tl. Since all terms on the right side of eq. (^) satisfy tl, it follows that t satisfies 
this property. 

To establish t2, we first show that the extensions Tfe[g] have an almost homogeneous 
scaling under rescalings of the metric in the sense of t2. To see this, we consider a 
term C ■ a*u in the expansion ([77|). By Theorem 4.1, the curvature term C will scale 
as under a rescaling of the metric by A^. On the other hand, for the term a*u, 
since is just the inverse of the exponential map at x, a rescaling of the metric will 
correspond precisely to a coordinate rescaling by a factor of A in the distributions u. By 
Lemma ^A\ , these distributions scale like A^""^ up to logarithmic corrections under such a 
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coordinate rescaling. Therefore, each individual term in formula ( [77| ) for has an almost 
homogeneous scaling with degree d under rescalings of the metric. On the other hand, 
the almost homogeneous scaling of under a rescaling of the metric can be proven by a 
argument similar to the proof that is local and covariant. Consequently, we see that t 
satisfies property t2. 

It also is relatively straightforward to prove that each occurring in eq. ([79D satisfies 
properties t3-t5. We know that is a finite sum of terms of the form C{x) ■ a*-u(?/), with 
C{x) a polynomial in the curvature and its derivatives. Since C{x) is smooth in x, we 
have 



WF{C ■ a*u) c{{x,Z[^yh;^i,[^]'h;...; 6 



kr, 



h-...- a,(e„), kn) G WF(m)|. (80) 

Here, we have written y = (^i, . . . ,^„) and each denotes a point in a convex normal 
neighborhood of x, and not the Riemannian normal coordinates of that point. (This makes 
a difference here, since we are considering variations in x.) In eq. (pOD, da^/d^i denotes 
the matrix of partial derivatives of with respect to C,i at fixed x, and da^/dx denotes 
the matrix of partial derivatives of ^^(Ci) with respect to x at fixed ^j. However, at = x 
we have dax{^i) / d^i = —dax{^i)/dx, since moving C,i infinitesimally away from = x has 
the same effect on ax{C,i) as moving x infinitesimally by the same amount in the opposite 
direction. It follows that if [x, ki, . . .]x, fc„+i) G WF(C ■ a*u), then Y^ki = 0. This means 
precisely that WF(C ■ a*u) fA„+i-L T{An+i), i.e., the microlocal spectral condition, t3, is 
satisfied. Similarly, by using the fact that C {x) is a polynomial in the curvature and is 
the inverse of the exponential map — so that both C{x) and have appropriate smooth 
and analytic dependence on the metric — together with the fact that u is independent of 
the metric, we find that the smoothness (t4) and analyticity (t5) conditions are satisfied 
by Tk. 

Thus, we would be done if our expression (^) for t corresponded to a suitably conver- 
gent power series. However, as already noted in remark (1) at the end of subsection 4.1, 
this is not the case, i.e., the remainder term, r^, in eq. (|79D is not expected to converge 
to zero in any sense useful for our purposes as m oo. Therefore, in order to prove 
that t satisfies properties t3-t5, it is necessary to explicitly analyze the remainder term 
Tm- This is technically quite cumbersome, since essentially the only thing useful that is 
known about is that it is the extension to A„+i defined by eq. (|78D of the expression 

given by eqs. (|56D and (|57|) . Equation (|78D expresses as a weak limit of distribu- 
tions whose wave front set properties are known, but wave front set properties are not 
preserved under weak convergence, so we must show that the sequence (|7^) converges in 
a suitably strong sense to enable us to prove that satisfies properties t3-t5. This will 
be accomplished in the proof of the following proposition, which — as will be explained in 
the remark following the statement of the proposition — will complete the proof that the 
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extensions t satisfy the properties t3-t5. 



Proposition 4.1. Let g^''-' be a smooth one-parameter family of smooth metrics, and 
let rm{s,xi, . . . ,Xn) denote the remainder term in eq. ([7DD, viewed as a distribution on 
M X f/". (Here m = d — 4n, where d is the scaling dimension of t.) Then the wave front 
set of Vm satisfies 

WF(r„) tMxA„± T(R X A„), (81) 

where the notation "_L" was introduced below eq. ([I6l). Similarly if g*^^^ is an analytic 
one-parameter family of analytic metrics, then (|81| ) holds for the analytic wave front set. 

Remark. If we choose g'-'^^ = g for all s, the above proposition implies that satis- 
fies the microlocal spectral condition t3. The proposition also implies that satisfies 
the smoothness and analyticity conditions, t4 and t5. In fact, the proposition asserts 
a somewhat stronger version of these conditions, as it shows that the wave front set of 
rm{s,Xi, . . . , Xn) not only cannot contain any points of the form (s, p; x, ki] . . . ; x, kn) with 
^ fcj 7^ but it also cannot contain any such points with p 7^ 0. Since each has already 
been shown above to satisfy t3-t5, it follows that t satisfies t3-t5 if Tm does. Thus, our 
construction of time ordered products satisfying properties T1-T9 of section 2 will be 
completed once we have completed the proof of this proposition. 

Proof. We will give the proof only for the analytic case; the proof for the smooth case 
is similar, though somewhat simpler because the estimates needed to establish the wave 
front set properties are simpler in nature in the smooth case. As before, we proceed by 
induction in the number, n, of variables (xi, . . . , Xn) on which depends. We inductively 
assume that the analytic wave front set version of eq. ( |8TD holds for all that depend on 
n or fewer variables. By a slight generalization of the proof given above that satisfies 
t3-t5, it can be shown that if g^*-* is an analytic one-parameter family of analytic metrics 



then each also satisfies eq. (|8TD . Consequently our inductive hypothesis implies that 



t{s,xi, ...,Xn)= t[g'^'*^](a;i, . . . ,Xn) satisfies 

WF^(t) rMxA„±T(Mx A„), (82) 

From the distributional coefficients t depending on n or fewer spacetime arguments 
and the real parameter s we obtain, by our inductive constructions, the distributional 
coefficients t°(s, xi, . . . , Xn+i) = t°[g^'*-'](xi, . . . , x^+i) depending on n + 1 spacetime argu- 
ments and the real parameter s. These distributions are defined everywhere in M x t/"^^, 
except for (M times) the total diagonal, A.„+i. Their analytic wave front set WF^(t°) is 
therefore a subset of T*{R x {U"'+^ \ A„+i)). By essentially the same arguments as given 
in [0, Sec. 7] (modulo a straightforward modification of those arguments with regard to 
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the additional parameter s on which the t° depend), the distribution t^{s,xi, . . . ,Xn+i) 
can be expressed as a finite sum of terms of the form 

t{s,{xi}i(zi)t{s,{xj}j(zic) Y\ Hpis^Xi^XjY'K (83) 



on each of the open sets Cj introduced in eq. (|TH]). Here, the aij are certain natural 
numbers, / is a nonempty proper subset of the set {1, . . . ,n + 1} and I'^ is its com- 
plement. (Note that since J is a nonempty proper subset, the expression (^) only in- 
volves the distributional coefficients t depending on n or fewer spacetime arguments.) 
Finally, Hf{s, xi, X2) = HF[g^^'>]{xi, X2) is the local Feynman parametrix introduced be- 
low eq. (|59D, for our analytic 1-parameter family of metrics. It can be seen by an explicit 
calculation that 

WFa{Hp) rMxA2^T(Mx A2). (84) 

For each of the sets / described above, let us define a projection map ttj from M x U'"''^'^ 
to M X t/'^' (with |/| the number of elements in /) by 

TT/ : (S, Xi, . . . , Xn+l) (S, {Xjjig/). (85) 



Using the rules for calculating the analytic wave front set of products of distributions [[1^] , 
we find from eq. (^) that the analytic wave front set of t° restricted to the open sets 
M X C/ is estimated by 

WF^(t°) txc.C (7r|WF^(t) U {0}) + (7r|. WF^(t) U {0}) 

+ E EKM}WF^(^^)U{0})cT*(Mx(f/"+^\A„+i))- (86) 

If we now take the closure in T*(]R x f/"+^) of the sets on both sides of the above relation 
(we denote this closure by an overbar), take the union over all /, and use eqs. (^^ and 
then we obtain 



WFA(tO) fMxA„+,^ T{R X A„+i). (87) 

From the properties of r^, it then follows that r^(s, Xi, . . . , x„+i) also satisfies the 
same condition, i.e.. 



WF^(r^) rMxA„+,±T(Rx A 



Note that eq. (^) imposes a nontrivial restriction (beyond what we already know) on the 
wave front set of r^. Our aim is to show that (^) continues to hold for the extension. 
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Tm- In order to simplify the discussion, we will show here only the weaker result that, for 
a fixed metric g, r^ixi, . . . , Xn+i) satisfies 

WF^(r„) U„+i^T(A„+i), (89) 

However, the arguments can be generalized straightforwardly to prove (^1]) in n + 1 vari- 
ables for an analytic one-parameter family of analytic metrics g*^^-*. 

As above, we choose relative coordinates {x, y) around the total diagonal. We will 
identify the point x E M with its coordinates in some chart, and we will identify y = 
(^1; • • • 5^n) with its Riemannian normal coordinates relative to x, so that the diagonal 
corresponds to y = 0. With this choice of coordinates, we identify t° (and likewise, r^) 
with a distribution defined on X x (F \ 0), where X is an open set in and Y is an 
open neighborhood of the origin M^". Let xo be some fixed point in X. It is possible 
to construct a sequence of smooth functions of the form (pj^{x,y) = (f)'j^{x)cl)'l^{y), where 
(f)'j^ G C^{K') is 1 in a neighborhood of Xq, such that ip"^ vanishes in a neighborhood of 
and is 1 outside some larger neighborhood K" , and where ipN satisfies the estimate 

IS^+'^^^I < Clfl+i(X + 1)1^1 V|/?| <X = 1,2,.... (90) 

If / is a test function with support sufficiently close to (xo,0), then the extension is 
defined by eq. ([7^). For our purposes, it is convenient to make the choice tD^'-''> = (0Ar)2^, 
where the subscript 2^ means the pull-back by the function {x,y) —>■ {x,2^y). 

In order to show WF A^rm) -L T{An+i), we must demonstrate that {xq, ko,0,pQ) is in 
the complement of WF A{rm) whenever ko ^ 0. It is not difficult to see that this will 
follow if we can show that it is possible to choose K = K' x K" so small that 

\ie^r0^ikp)\<2-^/'C''^\iN + l)/i\k\ + \p\)f (91) 

for all {k,p) in some conic neighborhood F of of (/^cPo) and for all natural numbers N 
and j. Here, 9^ G C^{K) is the cutoff function defined by (f)]\i{x,y) — (f)]y{x,2y). Note 
that the support of 6* at does not intersect the submanifold X x {0}, and that the sequence 
of cutoff functions 9n is again bounded in £\K) and satisfies the estimate (PD|). 

In order to analyze the Fourier transform on the left side of (pi]), we observe that 

{e^Tl{k,p) = 2-'-W;d^,-,{k, 2-^p), (92) 

where (r^)2-j denotes the pull-back of the distribution by the map (x, y) — > (x, 2^^y), 
and where the factor 2"^"-' is due to the fact that transforms as a density. Recalling 
our choice m = d — An, we can write the quantity on the right side of this equation as 

2-^ In 2ye^fi2-\ k, 2-^p), (93) 
I 
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where the ipi G x X x (Y \0)) were defined in eq. (|56D, and where the Fourier 

transform is with respect to the variables x and y. 

We now claim that for any closed conic set F in M x M'* x M^" not containing elements 
of the form (0,0, p) there is a neighborhood Kq C M x X x F of (0,a;o,0) such that for 
all/ 

WF^(V^°) n {Ko X F) = 0. (94) 

To prove this, we decompose ipf into simpler pieces, whose analytic wave front set is 
either known by the induction process or can be determined by elementary means. For 
this, we shall define a family of analytic metrics depending analytically on parameters 
s = (f,/i, x) & Pi X P2 X = P, where Pi is a small neighborhood of 1 in M, P2 is a 
small neighborhood of in M and where P3 is a convex normal neighborhood in M with 
respect to g. In order to define this family, let Xx,tJ. be the diffeomorphism which shrinks 
the Riemannian normal coordinates ^ of a spacetime point about the point x G P3 by a 
factor of fi. In terms of this family of diffeomorphism, our family of metrics is given by 

g(^) = {vfxr'xl^^g. (95) 

This is a real analytic family of analytic metrics (but with s now ranging over the 6- 
dimensional parameter space, P, rather than over M). We have already established that 
the analytic wave front set of the distribution t° on P x (f/"+^ \ A,„+i) satisfies ( ^7\ j (with 
M replaced by P). In order to relate ipf to we let R^"^^ be the map from test functions 
on M to smooth functions on M given by 



(p(-)/)(A) = ^ [\l 
ml Jo 



~f,rf<^+'\Xf,)dfi. (96) 



Furthermore, we set 



D^^^ = ^{vd/dv + dy. (97) 

Note that if / is a smooth function on M with compact support, then we have supp(*P''™-'/) C 
supp(/), where ^P^™-) denotes the transpose of R^"^\ Thus ^P^*") has proper support. It 
now follows straightforwardly from the definition of ipi that we can rewrite the action of 
the distributions ipi on test functions / G C^(]R x (f/"+^ \ A.„+i)) as 

^?(/) = iff) [(*^«5( • - 1)) ® iCRt^^ ® K...x.Jf)] , (98) 

where the subscripts on the operators indicate on which of the variables (f , /i, xi, . . . , Xn+i) 
they act, and where j*t^ denotes the pull back of by the analytic map 

J : {v, /i, Xi, . . . , x„+i) (f , /i, X = Xi, Xi, . . . , x„+i) G P X f/"+\ (99) 
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The analytic wave front set of ipf can now be estimated from eq. (^) using our knowledge 
of the analytic wave front set of eq. (^Tf), together with the rules for calculating the 
wave front set of a distribution under composition with distribution kernels [|T^, Thm. 
8.5.5], and under pull-back by analytic maps 0, Thm. 8.5.1]. For this, we only need 
to know the following additional facts: (i) The action of an analytic partial differential 
operator, such as D^''\ does not enlarge the analytic wave front set and (ii) the analytic 
wave front set of the distribution kernel of R^"^^ (viewed as a bidistribution on R x M) 
does not contain any elements of the form (A, 0;/i, p). The first statement is proven in 



1^, Thm. 8.4.7], and the second statement can be checked directly. 
This information suffices to conclude from eq. (PSI) that 



WF^(V'?) rMxA„+,^ T(M X A„+i). (100) 

If R X (f/"+^ \ A„+i) is identified with subset of M x X x (F \ 0) via the above choice of 
coordinates, then this means that 



WF^(V>0) tMxXxW^ T(M X X X {0}). (101) 



As a consequence, the open set T*(]R x X x F) \ WF^('?/'°) contains a set of the form Ko x F 
as claimed in eq. (|9^), provided that Kq is chosen to be sufficiently sharply concentrated 
about the point (0,Xo,0). 

We now blow up the sequence of cutoff functions Ojy G C^{K) to a bounded sequence 
cutoff functions in C^{Ko) which still satisfy the inequality (0), and which we shall 
denote by the same symbol. It then follows that with these cutoff functions, 

\e^^{p,k,p)\ < c^+i((x + i)/(|p| + \k\ + \p\)f (102) 

for all {p,k,p) G F, provided the support Kq of 9^ is sufficiently sharply concentrated 
near (0,Xo,0), where the tilde denotes the Fourier transform in R x R"^ x M^". With our 
new choice for 9^, we can write ( pSD as 

{0^rOm{k,p) = {27r)-'/'2-^J2(^\n2y [ 9^^{p,k,2-^p)e~''-''' dp, (103) 

where the sum is finite. Now the cone F can be chosen such that (p, k, 2~^p) G F for all 
points {k,p) in the cone F, all p and all j. Thus we can use ( |102|) to estimate 

\{9^Tl,{k,p)\ < 2-^/'C'^iN/\k\r~\ (104) 

For {k,p) in the cone F it holds that \k\ > e\p\ for some e > 0. This enables us to estimate 
the above expression further by 

< 2-^^^C^{N/{\k\ + \p\)f-^ (105) 

for all {k,p) in F and all natural numbers N and j. This is what we wanted to show. □ 
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Remark. The distributions t have now been shown to have an analytic dependence 
on the metric in the sense of condition t5. This makes it possible to establish certain 
analyticity properties of the distributions u in the scaling expansion for the t, as we will 
now show. We know that if (s, p; Xi, fci; . . . ; Xn+i, fcn+i) is an element in WF^(t), then the 
element (xi, ki, . . .; x„+i, A;„+i) must necessarily be in the set given by eq. (H). By 
our scaling expansion, we know that the distributions u are given in terms of s derivatives 
of t (evaluated at s = 0), so we can calculate WF a_{u) from WFyi(t) by the rules for the 
analytic wave front set under restriction and differentiation. It is straightforward to see 
that this gives 

WF a{u) C |(^i, ki]...] ^n, K) e T*M^'' \ {0} I 3 decorated graph G{p) in (M^ -q) with 
vertices 0, ^1, ^„ such that /cj = p^— Pe Vz|, (106) 

e:s{e)=i e:t{e)=i 

where we use the graph-theoretical notation introduced in T3, and where r] is the Minkowski 
metric. 



5 Conclusions and Outlook 

In this paper, we have given a construction of local, covariant time ordered products of 
an arbitrary number of local Wick powers. These local time ordered products were shown 
to satisfy properties T1-T9 of section 2. They therefore fulfill the assumptions of the 



uniqueness theorem of our previous paper |jT2|, Thm. 5.2]. Consequently, for any given 



polynomial order in the free field, any other prescription for defining local time ordered 
products with the same properties will differ from the prescription given in the present 
paper by products of local curvature terms and lower order time ordered products, as 
specified precisely in our uniqueness theorem. Although in this paper we considered only 
a massless Klein-Gordon scalar field, our results can be generalized straightforwardly to 
allow mass, and we do not anticipate any difficulties in generalizing our results to fields 
with higher spin. Largely for notational simplicity, we also restricted consideration to time 
ordered products of Wick powers that do not contain derivatives of the field, but it should 
be straightforward to generalize our construction to allow Wick powers of differentiated 
fields (subject only to the caveat of footnote |^). 

An important tool in our analysis was the scaling expansion introduced in subsection 
4.1 for the distributions t appearing in the local Wick expansion of time ordered products. 
In essence, this scaling expansion gives corrections to the "scaling limit" of |Tl[], organized 
in powers of the curvature (and the dimensionful parameters, if any are present). The 
scaling expansion generalizes to arbitrary t in the local Wick expansion the usual "short 
distance expansion" for the Feynman propagator (see remark (3) at the end of subsection 
4.1). Although we restricted consideration here to the distributions t, a similar scaling 
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expansion will exist for any local, covariant field which satisfies appropriate wave front set 
and scaling properties. The properties of the general scaling expansion for local covariant 
fields is currently under investigation. 

The results of this paper essentially complete the analysis of the existence, unique- 
ness, and renormalizability of the perturbative expansion of nonlinear quantum fields 
(with polynomial self-interaction) in curved spacetime. It is natural to ask whether an 
"exact" formulation of nonlinear quantum field theory in curved spacetime can be given. 
The Wightman axioms and other similar systems cannot be straightforwardly generalized 
to curved spacetime on account of their essential usage of Poincare invariance and the 
existence of a preferred vacuum state. We are currently investigating the possibility that 
the notion of a local, covariant quantum field (together with suitable microlocal spectral 
conditions, etc.) may enable one to give a useful formulation of axiomatic quantum field 
theory in curved spacetime. 
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A Smooth and analytic variation of distributions 

A key requirement that we impose on our definition of Wick polynomials and their time 
ordered products is that they have appropriate smooth/analytic dependence on the space- 
time metric. The purpose of this appendix is to elucidate the notion of smooth and 
analytic variation of distributions. 

To begin, let A be a manifold and for each s G M let m*-*-* : A — > C be a smooth (i.e., 
C°°) function. It is useful to view u^^^ as a map -u : M x A C. We say that u^'^^ varies 
smoothly with s if the map u is smooth. Note that this requirement of (joint) smoothness 
of the map u is stronger than the possible alternative requirement that u^'^\x) be a smooth 
function of s for each fixed x E X. This latter notion of (separate) smoothness in s would 
not be a natural one in the context of this paper for the following reason: We consider 
one-parameter families of spacetimes (M, g*^*^) and there is no natural way of identifying 
spacetimes with different values of s. However, the notion of separate smoothness is not 
invariant under diffeomorphisms ip^'^^ : A ^ A that are (jointly) smooth in (s,x). 

Now, for each s e M let m^*^ G P'(A), i.e., u^^^ is a distribution on A. We wish to 
define a notion of smooth variation of m*-'^-' with s that corresponds to the notion of (joint) 
smoothness of functions as defined in the previous paragraph. To do so, it is useful to 
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view u^^^ as a distribution, m, on M x X. The basic idea of our definition is to require u to 
be "not any more singular than each m*^'^^". One possible way of implementing this notion 
would be to demand that the wave front set of u be contained in the wave front set of u^^^ in 
the sense that WF(u) C {(s, p; x,k) | p = 0, {x, k) E WF(m'^''^)}. However, this definition 
is unsatisfactory for the following two independent reasons. First, the requirement that 
p = is too strong in that it would, in particular, require the singularities of u^^^ to 
"remain in a fixed location in X" as s is varied. This would not be invariant under a one 
parameter family of diffeomorphisms 'ip^'^^ : X ^ X that are (jointly) smooth in (s,x). It 
should be noted that the distributions, u^'^\ of interest in this paper have singularities on 
the light cones of g*^*-* and, hence, their singularities cannot "remain in a fixed location" for 
non-conformal variations of g. Consequently, we shall not require p = in our definition. 
Second, if m*^^-* happens to be "less singular than normal" for some value of s, then under 
the above proposed definition, u would fail to vary smoothly with s even if, in a naive 
sense, its variation with s was perfectly smooth. For example, u^'^\x) = s5{x) would fail to 
be smooth at s = because WF{u^^^) = but WF(m) includes points with s = 0. For this 
reason, we will define a more general notion of smoothness with respect to an arbitrary 
specified family of cones C^'\ (Here, a cone C is a subset of T*X\{0} having the property 
that if {x, k) e C, then (x, Xk) G C for all A > 0.) For the definition to be nontrivial, we 
must choose C^^^ so that WF(m('^)) C C^^^^ but we need not choose C^"^^ = WF(m(*)). 

Definition A.l. Let m*^*^ be a one-parameter family of distributions on a manifold X 
and let C^^^ be a family of cones. We say that u^^^ varies smoothly with s with respect to 
C*^^-* if the wave front set of the corresponding distribution m on M x X satisfies 

WF{u) C {(s, p; X, k) G T*(M X X) \ {0} | (x, k) G C^^^} (107) 

Remarks. (1) To illustrate the meaning of the above definition, let us consider the two 
extreme case, namely (a) when the cones are trivial, C^^^ = 0, and (b) when the cones are 
maximal, C^'^ = T*X \ {0}. In the first case (a) we immediately get that WF{u) = 0, so 
u is smooth jointly in {s, x). In the second case (b), it might appear that our smoothness 
condition is in fact empty. However, this is not the case, since eq. ( p.07| ) implies that no 
element of the form (s, p; x, 0) can be in WF{u). Thus, for example, if m = f ® with v a 
distribution on M and (p a distribution on X, not depending on s, eq. ( p.07| ) requires that 
V is smooth. 

(2) Let u^i^ and u^2'' be two families of distributions which are smooth with respect to 
cones C[^'' respectively €2'^^ Then the rules for calculating the wave front set of a sum 
of two distributions gives that the family u[^^ + ^2'^^ is smooth with respect to the cones 
Ci^UC2 \ Likewise, if {0} ^ C[^^ + C2 \ for each s, then the product u^i''u2^ can be defined 
for each s and defines a distribution jointly in (s,x). Moreover, the rules for calculating 
the wave front set of the product of two distributions gives that product family is smooth 
with respect to the cones C[^^ + €2^ . 



37 



The above definition allows us to define the notion of the smooth variation of a one 
parameter family, iu^^\ of continuous states on the algebras W(M, g^^^) of the spacetimes 
(M, g^*^): We say that cj'^*^ varies smoothly with s if each of the n-point functions w^*-*^ of 
(^(«) — viewed as a distribution on M" — varies smoothly with s in the sense of Definition 
AA\ with C(") = WF(wi'^). A one-parameter family of fields : P(M") >V(M, g(")) 
in n variables will then be said to vary smoothly with s with respect to the cones C^'' C 
T*M"\{0}, if the corresponding distributions vary smoothly with s with respect 

to C^^ for all smooth one-parameter families of states uj^^\ Since continuous states on 
W(M, g) are precisely the Hadamard states whose truncated n-point functions are smooth 
for n 7^ 2 [jl3|, it follows that ^^'^^ will be smooth if and only if is smooth for 



all smoothly varying families of Hadamard states c^;*^*-' with smooth truncated n-point 
functions (n 7^ 2). This is the requirement that we have adopted in condition T4 of 
section 2. 

In [O, a different notion of "continuous variation" of ^^^^ was introduced in the 



case where ^^^^ is local and covariant. Our uniqueness theorems for Wick polynomials 
and their time ordered products used the hypothesis that they vary continuously under 
smooth variations of the metric. It can be shown that our requirement of smooth variation 
introduced here implies continuous variation in the sense of |]T^, so the uniqueness results 
of |12[ continue to hold under this replacement (as can also be verified more straightfor- 



wardly by simply repeating the proofs with the new hypothesis). It also can be shown 



that the construction of local, covariant Wick polynomials given in [|T2| satisfies our new 
smoothness requirement with C^^ = 

In order to define the notion of analytic variation of a one-parameter family of distri- 
butions, u^^\ on a real analytic manifold, M, we first recall the definition of the analytic 
wave front set. To begin, let m be a function on M™' which is a real analytic in a neigh- 
borhood of a point Xq in W^. Then it follows from Cauchy's integral formula, or rather 
its generalization to C™', that 

< + 1)1"! Va (108) 

in a neighborhood of Xq, where C is some constant and multi- index notation has been 
used. Conversely, if the above estimate holds for a function m in a neighborhood of xq, 
then u is real analytic in that neighborhood. 

Condition (|108|) can be formulated equivalently in terms of Fourier transforms. Namely, 



one can show that an estimate of the form ( |108|) holds if and only if there is a sequence 
Mat of compactly supported distributions equal to u in some open ball around Xq, which 
is bounded in the space £'{W^) of distributions of compact support, and which satisfies 



\uN{k)\<C'^+\{N + l)/\k\Y en. (109) 

This motivates the following definition. Let m be a distribution on X C W^. The analytic 
wave front set WFa(m) is defined to be the complement of the set of all points (xq, k^) in 
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X X (M™" \ 0) such that there is an open neighborhood U of xo, a conic neighborhood F 
of ko and a bounded sequence mat G S'{X) which is equal to n on f/ and satisfies ( |109| ) 
whenever G F. 

It is clear from the definition that u is given by a real analytic function in the neigh- 
borhood of points Xq such that WFa{u) contains no element of the form (xo,fco)- If 
/ : X — i> y is an analytic one-to-one map, then the analytic wave front set of the pull- 
back f*u is given by {{x,df^{x)k) \ {f{x),k) G WF^(-u)}. This makes it possible, via 
localization in analytic charts, to define in an invariant way the analytic wave front set of 
a distribution on a real analytic manifold X. 

In practice it is useful that u^- can always be obtained as a product of u and suitable 



cutoff functions, see |jT^ Lem. 8.4.4]: Let WF^(m) fl {K x F) = for some compact subset 
K of X and some closed cone F, and let xn be a sequence of cutoff functions in C^{K) 
such that for all a 

< Cjf + y\P\< N = 1,2,.... (110) 



Then un = Xnu is bounded in £'{X) and satisfies ( |1Q9| ) for all k E F. 

A one-parameter family of distributions, u^'^\ on a real analytic manifold, X, will be 
said to vary analytically with s with respect to the cones C(^) if eq. ( [T07|) holds with WF 



replaced by WF^i everywhere in that equation. The notions of analytic variations of states 
and fields can then be defined in complete parallel with the definition of smooth variation 



given above. This agrees with the notions previously introduced in [12|. 



B Properties of the distributions u in the scahng ex- 
pansion 

In the following proposition, we list some general properties which hold for any almost 
homogeneous distribution on M"*. The distributions u in our scaling expansion are par- 
ticular examples of such distributions, and therefore the proposition applies to them. In 
particular, combining the upper bound ( |1(J6 ) on the analytic wave front set with item (iii) 



in the proposition below, one can obtain detailed information about the analytic wave 
front set of the Fourier transforms u of the distributions u in the scaling expansion. This 
information suffices to establish that u is in fact an analytic function in a large portion of 
momentum-space, and that it is given by the boundary value of an analytic function for 
almost all momentum configurations (see |T^, Thm. 8.4.15] for the appropriate criteria 
when a distribution can be written as the boundary value of an analytic function). 

Proposition B.l. Let u G V{W^) be an almost homogeneous distribution of degree p, 
i.e., S^u = for some N eN, where = (^ y'd/dy' + p)^ . Then 

(i) WFa{u) C {{y, k) G T*M™ \ {0} I E y% = 0}. 
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(ii) u can be extended to test functions in Schwartz space and thereby defines a tempered 
distribution. 



(iii) u is again an almost homogeneous distribution, with degree m — p. Furthermore, 
we have 

e WF^(m) ^ {k,-x) eWFAiu) ifx^O,k=^0, 
X e supp(m) (0, -x) e WFa(u) if x 7^ O, 
k e supp(?2) ^ (0, k) e WFa{u) ifk^O. (Ill) 



Proof. Since u = 0, and since Sp has analytic coefficients, we have by |T^, Thm. 8.6.1] 
that 

WF^(n) C Char(5f ) = {{y, k) e T*M™ \ {0} | J^y^fc, = 0}, (112) 

where Char(P) is the characteristic set of a differential operator P, defined as the set of 
all {y, k) G T*R"' \ {0} such that p(?/, k) = 0, where p is the principal symbol of P. This 
proves (i). 

Let x+ S'lid X- be smooth functions on R with the property that x+ + X- = I5 cind 
such that X-i'f') = fo^' ^ ^ ^0 and X-(r) = 1 for r > 2ro for some Tq > 0. We can 
therefore write 

u{y) = x+{\y\)u{y) + x-{\y\My). (113) 

The first distribution on the right side of this equation is by definition of compact support, 
and therefore trivially a tempered distribution. Thus (ii) will follow if we can show that 
also the second distribution on the right side is tempered. In order to prove this, we 
first show that it is possible to write X-"^ "polar coordinates". For this, we note that 
WF(m) n A^*(§™-i) = by (i), which implies by ^ Thm. 8.2.4] that u has a well 



defined pull-back, v, to the unit sphere, S*" ^, in M™. It follows from this, and the almost 
homogeneous scaling of u that it is possible to write 



^ 1 poo n 

w(X-/) = V / / c.X-iry-P+'^-Hn^ rv{y)f{ry)drdii{y), (114) 



where (r, y) denote polar coordinates in M"^, dp is the standard measure on and the 
Cj are unspecified complex constants. Since w is a distribution on S™"^, there must exist 
differential operators Pi, . . . ,Pk on S™'"^ such that 

|t^(/^)| <C5^ sup \PMy)\ (115) 



Kk 
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for all test functions h G P(S'"^^). Moreover T^_(^r)r~^'^''"~^ In-* r is a smooth function on 
M which grows polynomially together with all its derivatives at infinity, and therefore is 
a tempered distribution. Combining these facts with eq. (|114D, we easily get the estimate 



\u{x-f)\<C sup (116) 

for some a, 6 G N and all / in Schwartz space, thus showing that X-'u is tempered. 

We come to the proof of (iii). That the Fourier transform of u scales almost ho- 
mogeneously with degree m — p follows directly from our definition. For the case of a 
distribution, u that scales exactly homogeneously with degree p, the remaining three rela- 
tions in eq. ( [111| ) correspond precisely to [jl4|, Thm. 8.4.18]. The proof given in |T^ , Thm. 



8.4.18] of the second and third relations in ( |111| ) can be applied without modification to 



distributions with almost homogeneous scaling. We therefore only have to prove the first 
relation in ( p.ll| ). 



Since the Fourier transform u is again a tempered distribution which scales homoge- 
neously up to logarithmic terms, the problem is symmetric and it therefore suffices to 
show that 

[yo, h) i WFA(n) ^ (fco, -yo) i a{u) (117) 

if I/O 7^ 0, fco 7^ 0. Choose compact neighborhoods K and K in \ of ?/o and such 
that 

WF^(u) n (is: X ^) = 0, (118) 



and a sequence of cutoff functions xn ^ C'^iK) such that (|110|) is valid for every a. 



We now estimate the Fourier transform oi Vf^ = Xnu in a conic neighborhood of y^. By 
Fourier's inversion formula and the convolution theorem, we have 

VN{-Xy) = [ uix)xNix - Xy) (Tx. (119) 



We now estimate expression ( |119| ) for \y — y^l < r and arbitrary A. For this, we consider 
two cases, first < A < 1, and second A > 1. We begin with the first case. Since m is a 
tempered distribution, we can estimate 



\VN{-Xy)\ <C y2 \^"9^XNix - Xy)\. (120) 

H<a.,\l3\<b 

Using (|110D, it is not difficult to estimate 



\y 



'd('xN{y)\<C{N + l), m,\a\<a,\(3\<b (121) 
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From this one obtains the estimate 



\vN{-^y)\ < C{N + 1) < CZ''''^\{N -M+ 1)/Xy'-"' (122) 

for all M < and all < A < 1. 

In order to estimate ( |119D also for A > 1, we use that u scales almost homogeneous up 
to logarithms. This enables us to write 

VM{-Xy) = A^+™ E ^ / n,{x)xN{X{x - y)) (Tx, (123) 

where the sum is finite and where Uj = Sj,u, with Sp = J^V^^/^y^ + P- Since Sp is a 
partial differential operator with analytic coefficients, we conclude that WF A{uj) is not 
bigger than WF^(m) and hence also has no intersection with K x K. We are now in a 
position to use exactly the same arguments as in the proof of fl^, Thm. 8.4.18] (modulo 
a trivial additional estimate due to the logarithms), to show that there holds the estimate 

\VN{-Xy)\ < C^-'^+^((iV - M + 1)/A)^-^^ ViV, A > 1, |y - yo\ < r, (124) 

for some natural number M. Together with (|123|) this shows that we have v^iy) < 
C^-^+^{{N - M + l)/\y\)^-^'^ for all y in the conic neighborhood 

{-AyGM'" I \y-yo\<r,X>0} (125) 

of — I/O and for some fixed M. This proves the proposition. □ 
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